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Abstract 

We derive the connected correlation functions for eigenvalues of large Her- 
mitian random matrices with independently distributed elements using both 
a diagrammatic and a renormalization group (RG) inspired approach. With 
the diagrammatic method we obtain a general form for the one, two and 
three-point connected Green function for this class of ensembles when matrix 
elements are identically distributed, and then discuss the derivation of higher 
order functions by the same approach. Using the RG approach we re-derive 
the one and two-point Green functions and show they are unchanged by choos- 
ing certain ensembles with non-identically distributed elements. Throughout, 
we compare the Green functions we obtain to those from the class of ensembles 
with unitary invariant distributions and discuss universality in both ensemble 
classes. 
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I. INTRODUCTION 



Of recent interest in random matrix theory has been the discovery and characterization 
of new universal behavior in the densities of, and correlations between the eigenvalues of 
various ensembles which may be relevant to the study of mesoscopic systems, disordered 
metals, and random surfaces PJ-fl. 

In this paper we present an investigation for the existence of universal behavior in the 
correlation functions of a particular class of random matrix ensembles, those comprised of 
real symmetric or Hermitian matrices with individually distributed random elements. A 
notable example of an ensemble in this class is the first random matrix ensemble considered 
by Wigner ||, i.e. real symmetric matrices whose elements assume the values ±a/VN with 
equal probability. Here we will refer to this class of theories as the "Wigner class" . 

It has been known for quite some time that in Wigner class ensembles the eigenvalue 
density is universal. More specifically, the averaged density of eigenvalues for an ensemble 
of matrices with individually identically distributed elements is the same for any chosen 
distribution, up to a scaling of the width of the spectrum. We define the averaged density 
of eigenvalues as 



where ip is an N x N matrix and (.) denotes averaging over an ensemble of such matrices. 
In general, the ensemble average for an operator, O, is obtained via the integral 



where P is a density function, to be defined in more detail later, which takes as its argument 
the matrix (p. For any Wigner Class ensemble one obtains the well known result 







p(A) = 

ttgt 




which is known as Wigner's "semi-circle" law. 
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In this paper we will compare our results for "Wigner class" ensembles to those of 
another class, ensembles specified by a matrix density invariant under a symmetry group 
transformation. These "invariant class" (or "trace class" in the terminology of Refs. 
ensembles have been shown to have non-universal eigenvalue densities. In particular, one 
finds the density of eigenvalues for an invariant class ensemble may be written 

p(A) = A P (A)V^A5 (1.4) 

with p(X) a polynomial which depends on the ensemble distribution in a complicated way 
||10|| . However, the two-point and higher, connected correlation functions for this class have 
been found to be universal in form in the limit that iV — > oo, where TV denotes the dimension 
of the matrices |l|-44|,|6]] . Their only dependence on the averaging ensemble is through the 
endpoint of the eigenvalue spectrum. 

In a previous work ||, Brezin and Zee suggested that the two-point and higher correlation 
functions for the Wigner class (in the N —>■ oo limit) were universal as well, having the same 
form as those for the invariant class. Here we will show that this suggestion was incorrect, 
and that the two-point correlation function is somewhat less universal than its invariant class 
analog, as it depends on two parameters of the distribution (to be defined below) rather than 
just one. In general the averaged fc-point connected correlation for a Wigner class ensemble 
will explicitly depend on the ensemble's first 2k moments, and will be a quantity of order 
Af-3fc/2+i to J\f~ 2k ~ 2 depending on which moments are present. 

Amusingly, the occurrence of universality is reversed between these two classes. While the 
eigenvalue density for the Wigner class is universal, the higher order connected correlation 
functions increasingly depend on the particulars of the ensemble distribution chosen. In 
contrast, in the invariant class the eigenvalue density is not universal, but the higher order 
connected correlation functions are. 

For the two-point connected correlation function in the Wigner class (which we calculate 
for the Hermitian matrix case in the following sections) we find an expression dependent on 
only the second and fourth moments of the averaging ensemble. Specifically, we obtain 
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where (3 = 1, or 2 depending on if one considers real orthogonal or Hermitian matri- 
ces, a 2 is the second moment of the ensemble distribution, r 4 the fourth moment, and 
±2cr are the endpoints of the eigenvalue spectrum, all to be more explicitly defined later. 
Within the Wigner class one may consider, for example, such arbitrary and disparate en- 
semble distributions as P^CVy) cx Q(c 2 /N — |y?r,-| 2 ), Pjy(<Pij) oc exp(— iV YX=i 9k\Vij\ 2k )^ or 
P N ((fij) oc l/2[A5((pij - c/y/N) + B5(<pij + c/y/N)} always obtaining Eq.flOD. 

If we set r 4 = in Eq. (|1.5|) we obtain the result for a Gaussian ensemble distribution. 
This is identical to the universal form of the two-point connected correlation that was ob- 
tained by Brezin and Zee for any distribution in the invariant class. The Gaussian case, of 
course, belongs to both the Wigner and invariant classes. However, in general the two-point 
correlation for the Wigner and the invariant class differ by an additional term, namely the 
second term in Eq. (|1.5| ). 

In Fig. and Fig.|| we have plotted the functions (ji — u) 2 p c (fi, v)\p=\ and (// — 
j/) 2 p c (/i, u) | ( g =1)T 4^ with the arbitrary choice a 2 = 1, r 4 = —2, and N = 100 for two 
specific values of \i. We see that the second term is capable of changing the shape of the 
curves substantially. Thus, two Wigner systems with identical eigenvalue densities may have 
very different two-point connected eigenvalue correlations. 

In a numerical work Kobayakawa, et. al. [|llj found poor agreement between numerical 



calculations of two-point correlations for specific Wigner ensembles and the universal form 
for the invariant class. It appears that, with the inclusion of the second term, Eq.( |1.5| ) is 
in good agreement with the numerical results [IT , see Figures [3] and |] which correspond to 



(19) and (20) in Ref.]. 

Recently, we have learned that Pastur, Khorunzhy, and Khoruzhenko |12[ have also 
calculated correlation functions for the general Wigner class, but using a completely different 
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technique. We feel it worthwhile to present an independent study, since the techniques we 
use may be of interest, and the form of the resulting equations we obtain elucidating. Our 



result ( |1.5| ) differs from theirs (Eq.(14) in Ref. [|T^]) in the sign of the second term. 

This article is organized as follows. In Sec.[TI| we employ a diagrammatic approach to 
calculate the averaged one, two and three-point Green functions, or resolvents, for general 
Wigner class ensembles in the large N limit. We then outline how all higher point Green 



functions may similarly be calculated. In Sec.fTJ we demonstrate a completely different 
approach to the same problem, namely a renormalization group (RG) inspired procedure for 
calculating the same Wigner class Green functions. Naturally, the methods of both sections 
return identical results. 

Since from the diagrammatic standpoint this problem takes the form of a large N, zero 
dimensional field theory calculation, it is straightforward to determine the relevant diagrams 
and sum them. The RG method, on the other hand, takes more work and proves less 
transparent. However, the RG method makes more apparent the freedom one has in choosing 
model ensemble densities without changing the resulting correlation functions. 



II. CORRELATIONS VIA DIAGRAMS 

The use of diagrammatic methods in the study of random matrix theories is not new. 
Such techniques have been successfully applied in both invariant class and Wigner class 
calculations |p|,[l3| with applications to the latter being limited to a few specific but important 
ensembles, namely the ones with Gaussian distributed elements. 

Here we will show one how may calculate via diagrams the hierarchy of correlation 
functions for the Wigner class in general. The known results for Gaussian ensembles will 
be obtained as a special case. We will begin by explicitly calculating the one and two-point 
functions before considering the general n-point case. 

Let (p be an element of an ensemble of dimension N matrices with individually identically 
distributed elements. We limit our consideration to Hermitian matrices, however the method 
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may be extended to other sets such as symplectic or real symmetric. Also, we will initially 
take the case of an even distribution, and save treating the inclusion of odd moments until 

sec. urn 



We choose to define the moments of the matrix element distributions so that the eigen- 
values of any matrix will be of order iV -1 with mean zero and range of order N°. If we write 
<p = X+iY, with Xj = Xp and Yj = —Y? this is accomplished by defining: (Xj) = (Yj) = 0, 
((X]) 2 > = (1 + 5 l J )C 2 /2N, ((Yj) 2 } = (1 - 5})C 2 /2N } etc. so that (if)) = 0, (<p)<p*) = C 2 /N, 
(((p^ip 1 *) 2 ) - 2(C 2 /N) 2 = C4/N 2 , etc., where each C n is some constant of order N° and 
(.) denotes averaging over the ensemble. Note that in the case of a Gaussian ensemble all 
moments but C 2 equal zero. 

In general the kth moment is C^jX k l 2 . The inverse powers of N appear as a result 
of our requirement that the eigenvalue spectrum be finite in width as N becomes large. 
Alternatively we might have defined the ensemble density as having moments C 2 ,C^, . . .Ck, 
and the matrices of the ensemble as having elements (fj = X % -/\/N +iYJ/ \/N, with the same 
effect. We have opted, however, to define the distribution as above in order to facilitate our 
later counting of the powers oil/N associated with each diagram. 

We define the matrix 



G N (z)) = I — ] , (2.1) 



.1 



and the Green function, or resolvent, as 

G N (z) = ^TrG N (z)y (2.2) 

We are interested in G(z) = linijv^oc Gn(z) from which we may extract the averaged eigen- 
value density via p(/i) — ^ImG(/i + ie) . 

Similarly we define the two-point Green function as 

G N (w,z) = ^TtG n (z)^TtG n (w)Y (2.3) 

In the large N limit G^(w, z) — > G(w)G(z); thus we define the connected function as 



G cN (w, z) = G N (w, z) - G N {w)G N {z). (2.4) 

We see that this will be a quantity of order N~ 2 . The connected two-point correlation may 
then obtained via 

p c 0i, u) = ~[G e (+ t +) + G c (- -) - G e (+, -) - G e (-, +)], (2.5) 

where G c (±, ±) = G c (/i ± ie, z/ ± 2<5)| e 5 ^ . Three-point and higher Green functions may be 
defined, and their corresponding correlation functions obtained in a similar fashion. 

A. One point Green function 

To calculate G{z) we consider the untraced average of G{z) written as a power series 
expansion in 1/z, and then evaluate the ensemble average of each term in the series. Thus, 
we start with 



OC 



Gn{z)) = £ 



1 



n=0 Z 



Z 



(2.6) 



The average of any given term in the series may be obtained by summing over all the possible 
ways of taking the averages of groups of <^'s — pairs, triplets, etc. For example, the n = 4 
term of Eq. (|2.6|) is given by, 



1 // 1\ 4 \ 1,1.1 1 1 1.1 1.1 1 1.1 1 1.1 1.1.1.1.1 , . 
~ ( V 9- ) = -ip-(p-<p-<p- + -tp-<p-(p-(p- + -<f-(f-<p-(f- + -ip-ip-ip-ip-, (2.7) 

z \ \ z J z z z z z z z z z z z z z z z z z z z z 



where the constituents of a group average are indicated by having either a tilde or hat above 
them. Given our choice that (</?*■) = 0, single and triple ip averages do not contribute to this 
particular term. It should be emphasized that this procedure entails no approximation 

In order to evaluate ( |2.6|) by this procedure we consider it as a diagram expansion and 
use the double line formalism of 't Hooft fT4"| , where the lines in the diagrams correspond 
to the indices of the matrices being averaged. Figure || shows the diagrams associated with 
the terms in Eq. (|2.7Q . Each single line corresponds to a Sj/z and each double line connects 
a (fj to one or more others in an average. To borrow from the particle physics literature we 
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might call these "quark propagators" and "gluon propagators" respectively. Then two <^'s 
averaged together correspond to an emitted and absorbed gluon. The higher moments of 
a particular distribution correspond to gluon "interaction vertices" such as the four gluon 
vertex depicted in Fig. |](d) [P>1 - 

The great advantage of the double line formalism is that it allows one to count the powers 
of N associated with each diagram with ease. Each continuous closed line represents a sum 
over the values of an index, and contributes a factor of N. Thus, non-planar diagrams, such 
as Fig. |(b), are suppressed by powers of 1/N relative to equivalent planar ones since a 
non-planar diagram always has fewer closed index sums than a planar diagram having the 
same "interactions". 

The nature of the gluon-gluon interaction vertices affords another simplification in the 
large N limit. Each kth order gluon interaction vertex, such as the fourth order one in Fig. 
|5](c), contributes a factor of Ck/N k / 2 (in this way of counting the powers of 1/N associated 
with the gluon propagators are included). However, any such vertex has only two indices 
which may be summed over. Figure |6](a) depicts a four-gluon vertex with lines labeled by 
their index. If such a vertex appears in a diagram, that diagram will have a contribution of 
order N~^ k ' 2 " 1 ' , or less, to the total expansion. 

This is the crucial, and simple, reason why the one-point Green function for this class 
of ensemble is universal. In the large N limit the expansion of Eq. (|2.6|) only depends on 
diagrams containing averages of pairs of matrices. In particle physics language we would 
say that the interactions between the gluons becomes weaker and weaker as N approaches 
infinity. Thus for large N, regardless of the particular features of a distribution, the one- 
point Green function and averaged eigenvalue density depend on only one number: the 
ensemble density's second moment. 

This is in direct contrast with the result for invariant class ensembles. Invariant class 
one-point Green functions are known to depend on the ensemble for which they are calcu- 
lated in a complicated way [ rOfl . They are completely non- universal. From a diagrammatic 



standpoint, the key difference between Wigner and invariant ensembles is seen in the gluon 
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interaction vertices. The diagram expansion entails summing over all possible ways of taking 
the the averages of groups of matrices in an expression. In the case of an invariant class 
problem, enumerating the diagrams to be summed involves perturbatively expanding the 
ensemble density around the Gaussian case jj. The perturbations, which are identified with 
gluon interaction vertices in diagrams, are invariant under the same transformation as the 
averaging ensemble. Because of this, a kth order gluon interaction vertex will have k indices 
which are summed on rather than just the two that an analogous Wigner type vertex would 
have. Such vertices may not be neglected in the large N limit, and by their inclusion the 
one-point Green function becomes dependent on the distribution in detail. Diagrammatic 
representations of a simple four-point vertex are displayed in Figure |6| in order to contrast 
their index structures. 

Returning to our Wigner class calculation, the diagram expansion for G(z) l j looks like 
that in Fig 0(a), in the large N limit. By introducing £(z)j as the sum of irreducible 
diagrams — that is, diagrams which may not be cut on a single quark line to produce two 
complete and independent ones — we obtain the "integral" equation, 



which is depicted in Fig. 0(b). It is obvious from the graphs that £(/z)j = ^)^{ z )^ thus Eq. 
( p.8| ) can be written as a simple quadratic equation. Solving this for S, and using the fact 
that £(z) = cr 2 G(z) yields the Green function 

G(z) = ^(z-y/z 2 -4a 2 ) (2.9) 



This of course corresponds to the eigenvalue density Eq. ( |1.3| ), a long known result. Since 
it only depends on a 2 and not on any other details of the ensemble distribution, it holds 
for any ensemble in the Wigner class, and therefore is known as a universal result. In our 
further results, however, we will find an increasing dependence on the particulars of the 
ensemble distribution starting with the two-point connected Green function. 
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B. Two-point connected Green function 



We write G c (w,z) as an expansion in 1/z and 1/w with two derivatives taken out ex- 
plicitly. 



where the subscript c denotes the connected part. In diagram language this expression 
is a sum of graphs that consist of two different quark loops interacting via one or more 
gluons. Notice that since we have taken out the two derivatives there are an equal number 
of </?'s and propagators in each trace. The action of each derivative is to sum over all the 
possible insertions of a quark-quark vertex on it's respective quark loop. The evaluation of 
the diagrams will be made easier by having taken out these derivatives. 

Consider first the case where there are no gluon-gluon interactions, only quark-gluon 
ones (namely, the Gaussian case), and where contractions are made only between the two 
traces. In this case then m and n must be equal. If we draw the diagrams so that one loop 
is inside of another then the planarity of graphs is easy to determine. There are n ways of 
connecting the first gluon from one loop to another and then only one way of drawing the 
rest without crossing. Each connection corresponds to averaging together a (p in one trace 
with a ip in the other and results in a factor of a 2 /N . In the double line formalism we see 
each gluon connection creates a new closed index loop and thus contributes a factor of N . 
These graphs sum to 



Next consider higher order gluon interactions, still only treating inter-trace connections. 
Each gluon-gluon interaction vertex of order k, corresponding to the averaging together of 
k </?'s split between the traces, contributes a multiplicative factor of Ck/N k / 2 while adding 
at most one additional summed index loop. Since G c (w, z) is of order N~ 2 , in this case only 
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one such graph is allowed, the one with four gluons, two on each loop, interacting via a four 
point vertex (k — 4). This graph, depicted in Fig. |(a), averages to 

r^d_d_ /J_\ 2 
2 dw dz \wz) 

with r 4 = C4, the fourth moment of the ensemble distribution. 

We must next consider contractions within the same trace. These consist of what could 
be termed "vertex" and "self energy" corrections. By a vertex correction we refer to the 
inclusion of gluons which connect two quark lines in the same loop, on either side of a group 
of one or more quark-gluon vertices. We immediately see that vertex corrections may be 
ignored since such a gluon turns two closed index loops into one, suppressing the diagram 
in which it appears by at least N~ 2 . By self energy corrections we refer to the inclusion of 
gluons emitted and re-absorbed between vertices. The inclusion of self energy corrections 
to the quarks is instantly achieved through the "dressing" of the bare quark propagators 
1/w and 1/z by taking 1/w — > G(w) and \jz — > G(z) everywhere, excepting the partial 
derivatives. 

Thus, we obtain the complete two-point connected Green function for Wigner class en- 
sembles in general, 



d d 



dw dz 



N 2 G c (w, z) = -— — log 1 - a 2 G(w)G(z) - -G(w) 2 G(zy . (2.11) 



T 



2 



For the case of Gaussian distributed matrix elements r 4 = and Eq. (|2.11| ) reduces to the 
well known result for the Gaussian unitary ensemble (GUE) as we would expect. In general 
though, we find a dependence on two parameters, the second and fourth moments of the 
distribution. 

C. Three-point and higher Green functions 

Here we will discuss the application of the diagram method to the calculation of n- 
point connected Green functions. In the previous section we saw that, before the dressing 
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procedure, the expansion for the general two-point function differed from that for the GUE 
by only one diagram. This diagram, once dressed, may be viewed as a correction to the 
GUE result. As we will see, the higher order connected functions lend themselves to a 
similar interpretation, and may be written 

G c n(zi, z 2} . . . , z n ) = G^(zi, z 2} . . . , z n ) + (non-Gaussian terms), (2-12) 

where the superscript (0) indicates the GUE result. From here on we will not explicitly 
calculate the GUE Green functions as they are well known in the literature (see, for example, 
fl3fl or H). Our main interest will be on the non-Gaussian terms and the diagrams that 
comprise them. 

The relevant, undressed non-Gaussian diagrams are straightforward to enumerate. By 
our conventions an n-point connected Green function for the GUE, as well as the diagrams 
in its associated expansion, is of order ]\f- 2n+2 _ One creates the non-Gaussian diagrams for a 
general ensemble by considering gluon interaction vertices; these have a few qualities which 
greatly reduce the number of relevant diagrams in the large N limit. 

First, a kth order gluon interaction vertex (k > 2) contributes a multiplicative factor of 
Ck/N k l 2 while only having at most two free indices associated with it which may be summed 
on. Including such a vertex in a connected diagram makes the diagram order ]\[- n - k / 2 + 2 or 
less. This fact limits the moments that the n-point function can depend on to those up to 
Cm- As we have seen the 2-point function depends only on C 2 and C4, the 3-point function 
will depend on C 2 -,C^ and Cq, etc. This fact also limits the number of vertices which may 
appear in any one diagram. In the 2-point case only a single k = 4 vertex is allowed; the 
inclusion of another creates a diagram of order N~ A or less. In the three-point case we may 
have a single k = 6 vertex, or up to two k = 4 ones. 

Second, any two quark loops connected via a gluon interaction vertex may not have 
any other connections between them to leading order. The simplest case where this may 
be observed is in the single non-Gaussian diagram of the two-point function, depicted in 
Fig.[8](a). If we add a single gluon connecting the two loops we obtain Fig.|8|(c). This new 
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graph is now down by N^ 1 because we have added a gluon contributing a 2 /N yet created 
no new sum loops since the gluon interaction vertex requires the lines on either side of it 
to have the same index value. This effect persists in the general n-point connected Green 
function. 

As a demonstration of the use of these simplifying facts we calculate the three-point 
connected Green function to leading order. From our previous discussion we determine that 
the non-Gaussian diagrams of the three-point function will depend only on the C4 and Cq 
moments and the diagrams in their expansion will have either one 6-point, one four-point, 
or two four-point gluon interaction vertices [see Fig. 0(a)] . 

There is only one, undressed diagram containing a 6-point gluon interaction vertex which 
is relevant to order iV~ 4 . It corresponds to a term 



with v 6 = C 6 . 

In the case where we have a single four-point gluon interaction vertex we have two sub- 
cases which are not suppressed: one with the vertex connecting all three loops, and the other 
with the vertex connection two of the three. With the vertex connecting all three there is 
only one basic diagram (plus permutations), two gluons connecting to one quark loop and 
one each to the other two. This corresponds to, 



If the four-point vertex only connects two of the quark loops there are an infinite number 
of graphs. This occurs because, although once we have connected two of the three quark 
loops via a four-point gluon vertex there may be no other connections between them, either 
one or the other of the loops may have any number of non-interacting (Gaussian) gluon 
connections to the third. The resulting diagram expansion will be identical to the expansion 
of the two-point connected Green function summed over all possible ways to connect a third 
loop via a four-point gluon vertex. We use the distributive property of the derivative to 
immediately write this as 




(2.13) 




(2.14) 
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-h did2d3 u(i d2+ i d3 



log 1 



a 



r 4 1 



<-> 2) + (1 <-> 3U . (2.15) 



^2^3 / 2 Z 2 ^3. 

A note, relating to the interpretation of this problem as a field theory calculation, is that 
we would not be able to obtain this sum of diagrams so easily if it were not for the fact we 
are dealing with a static, zero dimensional field theory. In a field theory with "time" and 
"space" variables, quark lines would carry momentum which would be conserved at vertices. 
Then the derivative trick would no longer work. 

Finally, one can see that the inclusion of diagrams with two four-point gluon interactions 



is achieved by substituting G^i' for Gc23 Und i n Eq. fl2-15|) . Putting all the terms together, 



dressing them via — > G?j, and simplifying some, we find the three-point connected Green 
function is 



N 4 G cl23 = N 4 G^ 23 - d 1 d 2 d 3 ^r A G 1 G 2 G 3 {G 1 + G 2 + G 3 ) + ^G\G\G\ 
+ r 4 \G\ (G 2 d 2 + G 3 d 3 ) (log(l - a 2 G 2 G 3 ) - ^-G 2 G 2 3 ) + (1 <- 2) + (1 <- 3)] 1. (2.16) 



D. n-point connected Green functions for C n ^ 

Following the procedure above one may generate the entire hierarchy of connected Green 
functions. In general, the nth connected Green function will equal the GUE n-point function 
plus additional terms depending the ensemble distribution's moments up to C 2n and on 
powers of the one-point function and its derivatives. However, if the ensemble has an nth 
moment the form of the n-point Green function will be far simpler. 

Consider our previous example. In the three-point function calculation above we specified 
an even distribution for the ensemble. If we allow both even and odd moments of the 
ensemble distribution to exist we must also consider for this case diagrams containing three- 
point gluon interaction vertices. We define a third moment of the distribution: ((Pjif^iPj) = 
( W ( 4 — ^)$j)C3/N 3 / 2 (no sum on repeated indices). Then, three additional, undressed 
diagrams, depicted in Fig.|T0|(b), contribute to the three-point connected Green function. 
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The second two are order N~ A — the same as all the original terms in Eq. ( |2.16| ). The first, 
however, is order N~ 7 ^ 2 . In other words, this single diagram is order \Z~N larger than all the 
others. Thus, we have 

N 7/2 G cl23 = -ed x d 2 d z (G 1 G 2 G 3 ) + 0(iV- 1 / 2 ) (2.17) 

where c 3 = C 3 . 

For an ensemble with a non-zero nth moment (n > 2) one finds that the n-point connected 
Green function, to leading order, is 

n 

N 3n/2-i Gci ^ = Cn Y[(-di)Gi + 0(N- 1+1/2{n mod 2) ). (2.18) 
1=1 

Thus, depending on the particulars of the ensemble density, the leading order n-point con- 
nected correlation function can be a quantity of order JSf~ 3n / 2+1 through order jy -2 ™ -2 . This 
sort of dependence on the existence of various ensemble moments is in striking contrast to 
the case in Invariant class ensembles. For the Invariant class, 2-point and higher connected 
correlation functions depend on the ensemble density only through the endpoint of the eigen- 
value spectrum and the order in 1/N of the various quantities is, in general, independent of 
the ensemble density. 



III. RENORMALIZATION GROUP APPROACH 

In this section we demonstrate the use of a renormalization group inspired approach 
16| to calculate the connected two-point Green function for theories in the Wigner class. 



Further application to higher point functions follows in an obvious way. 

This method has been used previously, by Brezin and Zee 0, to derive the one-point 
green function for the Wigner class of theories, as well as by others to derive various quan- 
tities in invariant class theories [IT]]. In the interest of making the somewhat complicated 



derivation of the two-point function more transparent we will briefly review Brezin and Zee's 
calculation of the one-point function. 
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A. One-point Green function via RG 



The basic idea here is to write the (N+l) dimensional, matrix Green function in terms of 
dimension iV quantities and average over select elements. This generates a partial differential 
equation which may then be solved in the large N limit. 

To begin the derivation of G(z) we write the (N + 1) x (N + 1) Hermitian matrix <^v+i 
in the form 



( \ 



with (ppj an N x iV Hermitian matrix, v a complex N component vector, and \ a real scalar 
then the identity 

1 _ 1 



1 + (v\ 


{ 1 ) 


2 


V) 


z-X- 


(v\ 1 

\ 1 z—tp N 


V) 



Tr = Tr + V'"^/ (3.1) 



may be established. One can write this as a partial differential equation for G, 

TrG N+1 (z) = TrG N (z) + ^ \og(z - X - (v\G N (z)\v)). (3.2) 



We obtain a differential equation for G(z) by averaging Eq. ( |3.2| ) over the distribution 
governing <^n+i, and retaining terms of order iV° or larger. 
The term on the left hand side of Eq. ( |3.2|) averages to 

(N + l)G N+1 (z) = (N + l)G N {z) + (N + \)^G N {z) + ... 

Since Gn(z) is of order N°, to the desired order we may neglect all but the first term on the 
right. 

Next, we skip to the second term on the right of Eq. ( |3.2|) . Recall that the mean squared 
value of matrix elements is of order iV -1 . To order iV°, then, we may immediately set 
X = in this term. If we expand the logarithm we seen that we will be required to evaluate 
averages such as (v*VjVlvi . . .v*v q ). Consider the term with four t>'s 

(v*VjV* k vi) = A(5 i:j 5 k i + S u S kj ) + B5 ijkl 
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In the large N limit A terms will dominate, with A = (u*uiufu2) = (vlvi) 2 = (a 2 /N) 2 In 
general every average of a product of v's will be dominated by second moment terms, and 
in fact we may replace (v\l/(z — (pp{)\v) with (a 2 /N)Gn(z) within the logarithm. 

We must be more careful in averaging the first term on the right of Eq. ( |3.2| ) since we 
are averaging over the distribution appropriate for (Pn+i, and not tp^. Our experience with 
the second term, which we treated in the previous paragraph, indicates how free our choice 
of a matrix density may be. First, G(z) depends only on the second moment, a 2 , of the 
distribution. In fact we may set off diagonal elements to zero with some finite probability 
without changing G(z). Second, diagonal and off diagonal elements may even obey different 
distributions since all that is required by this procedure is that (x 2 ) — > as iV — > oo. 

Thus, if P((pij\ cr 2 /N, q 3 /N 3 ^ 2 , t 4 /N 2 , . . .), defines a particular distribution (with a 2 , q 3 , 
r 4 , . . . being its moments) then it is sufficient to use P{(fij', & 2 /N, 0, 0, . . .) in the calculation 
of G(z), that is, we may effectively set q,T,. . . equal to zero. After some straightforward 
manipulation we see that first term on the right of Eq. ( |3.2| ) averages to NG z; (l — -^j a' 
0(N- 1 ). 

Putting all this together, we find that in the limit iV — > oo, Eq. (|3.2|) averages to 



+ 



G(z) = -a 2 -^G(z) + A log^ - a 2 G(z)) (3.3) 

From dimensional analysis one can see that G(z; a 2 ) = (l/a 2 )G(z/a;l) and hence 
a 2 £zG{z) = -|(1 + z-§- z )G{z). Finally, we find that Eq. Q becomes 

G(z) = z^-G(z) + 2^- hg(z - a 2 G(z)). (3.4) 
dz dz 

This nonlinear equation has the unique solution 

1 



G{z) = —Az - V^ 2 -4a 2 ). (3.5) 



Naturally, this G(z) is identical to that in Eq. ( |2.9| ) and thus corresponds to the eigenvalue 
density ( |1.3| ). Now we will apply this method to the problem of deriving the connected 
two-point Green function. 
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B. Two-point connected Green function via RG 



As a starting point, we use Eqs. Q2.1J) and (|3.1|) to write 



Tr Gn+i(w)Tt Gn+i(z) = Tt Gn(z)Tt Gn(w) 

+ -^-TtG n (z) log(tw - X - (v\G N (w)\v)) 
d 

+—TiG N (w)\og(z-x - (v\G N (z)\v)) 
d 2 

+ dwdz l ° g ^ W ~ X ~ \ v \ GN y w )\ v ') lo s( z ~ X ~ (v\G N (z)\v)). (3.6) 

Now we must average over the distribution governing <Pn+i> keeping only connected terms 
up to order AT -1 . 

The term on the left of Eq. fl3~6|) is trivially averaged to 



(N + l) 2 G c{N+1) (w, z) = (N+ l) 2 G cN (w, z) + (N + l) 2 -^G cN (w, z) + . . . 

Since we are keeping only terms up to order iV -1 , and we know that G c n(w, z) itself is order 
iV~ 2 , then we may write this as 

N 2 G cN (w,z) + 0(N- 2 ) (3.7) 



Next we consider the second (and third) term on the r.h.s. of Eq. ( |3.6|) . Let us first 
obtain the connected average to leading order in 1/N. To leading order we may immediately 
set x — and take the average over v inside the logarithm where we see we may again replace 
{v\Gn(w)\v) with (a 2 /N)TrG(w). Thus, we have 

Now we complete the average over the remaining variables, here using the distribution 
appropriate to (p^ since corrections are of order iV -1 and still only concerning ourselves 
with leading order contributions. To leading order we need only average together TrG(z) 
and one of the TrG(w)'s in each series term. Then Gn(w) may be substituted for the 
remaining Tr(?(u>)'s. By summing the series, and utilizing the identity 



G N {w) 



1 

w — cr 2 Gjv(w 



0(iV- 1 ), 



we obtain 



d 

N—a 2 G cN (w,z)G N (w) 



(3.9) 



as the leading order average of the term in question. Note that this expression is already 
order iV -1 , therefore we need not consider any corrections to it as they will be order N~ 2 or 
smaller. 



On to the fourth term in Eq. ( |3.6| ). To the desired order we retain at most a single 
second moment in x, any higher moments or powers of moments may be neglected, yielding 



d 2 



dwdz 



00 1/1 ~ \ n /l 
£ —(-(v\G N (w)\v)) (-(v\G N (z)\v) 



\m,n=l 



ran \w 



+~ (fi + 1 ~<«|6»(«)l«>-) f 1 + t ~W6^M 



n=l 



n w' 



m=l 



m z" 



(3.10) 



In the first term of ( P-10 ) contributions to order 1 will come from one of two connected 
averages: either 



(^)(v k v^d(w)*G(z)t = ^Tr [G(w)G(z 



a 

N 2 

4 N 



or 



M^;) c GH}G(z)f = l-Y,G{w))G{: 

1=1 



(3.11) 
(3.12) 



where r 4 is the fourth moment of the ensemble distribution. Each term is of order iV l , 
thus either one or the other will occur once. The contribution from ( |3.11| ) is 



1 d 2 
N dwdz 



w — <j z G(w) z — a z G(z) \N 



1 d 2 



N dwdz 



G{w) - G{z) 



a i G{z)G{w) 



w — z 

1 d 2 a 4 G 2 (w)G 2 (z) 
Ndwdzl-a 2 G(w)G(z) 



(3.13) 



and the contribution from ( |3.12|) is 



_4 N 



1 d 2 

Ndwdz [w - a 2 G(w) z - a 2 G(z) \ N ^ 



J2G N (w)lG N (z)l 



-G 2 (w)G 2 (z) (3.14) 
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Treating the second term in (|3.10|) is far simpler since its average is unconnected and it 
is manifestly of order iV _1 . We may immediately average it to 



N zw 



1 - a iG{w) 



1 + 



°iG{z) 



a 



—G(w)G(z) 



(3.15) 



Finally, we come to the first term on the r.h.s. of Eq. (|3.6|) , we have a case where, as in 
the calculation of G(z), we must carefully average over the distribution appropriate to <^tv+i 
and not </?jv- Again we are led to consider how free our choice of probability distribution may 
be. Consideration of the previous terms has shown that G c n{w,z) will explicitly depend 
on only two numbers a 2 and r 4 , thus it is sufficient to average this term using the general 
distribution P(^; a 2 /N, t a /N 2 ). 



We find the term in question then averages to 



1-4, a' 



d 



+ 2t 4 



d 



N V da 2 dr 4 
From dimensional analysis one can see that 

1 



G nN (w,z) 



G cN (w,z; a 2 ,T 4 ) 



G cN (w/a,z/a; l,r /(a 



2\2\ 



Hence, the identity 



Using this we see that the first term on the r.h.s. of Eq. (|3.6|) finally averages to 



(3.16) 



Using expressions 
equation 



j|), (|3.14|) , (|3.15|) and (|3.16| ) we obtain the partial differential 



Y, d i (z i -2o 3 G i ) 

i=l,2 



N 2 G cl2 = 



a 2 dG 2 



+ r i G{G 2 2 )+0(l/N), (3.17) 



1-(t 2 GiG 2 

where we have used the obvious shorthand: G c u = G c (w, z), d\ = (d/dw), G± = G(w), etc. 
Note that the derivatives on the l.h.s. act on G c i 2 as well as on the parenthesized terms 
immediately to their right. 
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This rather complicated, non-linear equation is solvable. In general we may write its 
solution in the form 



N 2 G cl2 = d 1 d 2 F(G 1 + G 2 , G x - G 2 ) + ^(d) + h 2 {G 2 



(3.18) 



with F, hi, and h 2 unknown functions. We may immediately set hx t 2 = since G c i 2 must 
go like 1/zl 2 as Z\$, — > oo. We substitute this general solution into Eq. (|3.17|) and use the 
large N identityG/G" = 2a 2 G - z, where G' = (d/dz)G(z), 



d x d 2 



+ F(G 1 + G 2 ,G 1 -G 2 ) 



did 2 



2a 2 GiG 2 t A„ 2 „ 2 



+ t GfGi 



(3.19) 



l-2a 2 G 1 G 2 

The structure of the equation suggest that we make the variable substitution g± = 
-^(logGi ± \ogG 2 ). We equate the bracketed terms in Eq. ( |3.19| ) to obtain 

(3.20) 



^ 8 2n 2 p 9 + 
og + 1 — a 2 e9+ 



which is trivially integrable. Thus, 
N 2 G cl2 = -<9k9 2 



log(l - a 2 dG 2 ) - —G\G\ 



+ f(Gi/G 2 ) 



(3.21) 



where / is an unknown function. The known asymptotic behavior of G c i 2 requires / = 0, 
thus the unique solution to Eq. ( |3.17| ) is 



N 2 G c (w, *) = - A JL Lg [l - a 2 G{w)G{z) 



T 



G(w) 2 G(z) 



(3.22) 



By Eq. 



) we find the corresponding correlation function is 

4a 2 -fjiv + J (4a 2 -/i 2 )(4<x 2 - v 2 ) 



-1 d d 



log 



An 2 N 2 dfi du y 4cr 2 _ ^ _ y / (4a 2 -/i 2 )(4a 2 -z/ 



T [IV 

~2o^ 



4cr 2 - /x 2 a/4ct 2 - v 2 



(3.23) 

which corresponds, when differentiated, to the previously displayed Eq. (|01) . We simplify 
by a change of variables to sin^ = fi/2a, sin</> = vj2o which allows us to write 

-1 



Pc(M) 



16n 2 a 2 N 2 cos 6 cos 6 08 



d d ( l + cos(# + 
log- 



1 - cos(# - 



2r 4 

— - sin 26 sin I 

a 4 



(3.24) 
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We can compare this to the known result for invariant class Hermitian ensembles 

(f) -1 d d 1 + cos(fl + 0) 

Pc{ ,<P) 4nW W cos cos cfidOdct) 10g 1 - C os(# - 0) 1 } 

where a is the endpoint of the eigenvalue spectrum and depends in detail on the ensemble 
distribution. The Wigner class result is slightly "less universal" in that it has an additional 
term and depends on a second number, r 4 . As we have seen, higher order connected Wigner 
correlation functions will depend on more and more moments of their ensemble distribution. 



CONCLUSION 

For Wigner class ensembles, we have derived a general form for the one, and two-point 
connected correlation functions using both a diagrammatic and a renormalization group 
approach. Using the diagrammatic approach, we have also shown how one may construct 
higher order correlation functions, in the large N limit. By the RG inspired method we have 
demonstrated that results for the one and two-point functions are unchanged by considering 
some general cases where the elements of matrices in the ensemble are not identically dis- 
tributed. For instance, one may set elements to zero with some finite probability or allow on 
and off diagonal elements to have different distributions without changing the final result. 

In general the n-point function will depend on the ensemble's moments up to C 2n and 
be comprised of products and derivatives of the universal one-point Green function. We 
have found, though, in the case where a particular ensemble distribution has a non-zero nth 
moment the n-point connected Green function is particularly simple to leading order. This 
contrasts with the case for invariant class ensembles, where the one-point Green function 
depends on the ensemble in detail while higher order connected Green functions are universal 
in form. 
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FIGURES 

FIG. 1. The function (/x — z/) 2 p c (/i, ^)|^=o for the Wigner class with parameters f3 = 1, a 2 = 1, 
and iV = 100, plotted in the two cases r 4 = (solid line) and r 4 = —2 (dashed line). 

FIG. 2. The function (/i — z/) 2 p c (/i, i^)L=i for the Wigner class with parameters (3 = 1, a 2 = 1, 
and iV = 100, plotted in the two cases r 4 = (solid line) and r 4 = —2 (dashed line). 



FIG. 3. Fig. (19) from Ref. [11|. Dashed line is a numerically obtained smoothed connected 
correlation (/i — ^) 2 p c (A i > l/ )U=o f° r a rea l symmetric Wigner class ensemble with matrix elements 
chosen to be ±l/\/iV with equal probability, and N = 100. Solid line is (/i — z^) 2 times the 
universal connected correlation for real symmetric invariant class ensembles with N = 100, and 
spectral endpoints ±2, at ji = 0. 



FIG. 4. Fig. (20) from Ref. [11|. Dashed line is a numerically obtained smoothed connected 



correlation (fi — v) 2 p c {ii, v)\^ = i for a real symmetric Wigner class ensemble with matrix elements 
chosen to be ±l/y/N with equal probability, and iV = 100. Solid line is (fi — v) 2 times the universal 
connected correlation for real symmetric invariant class ensemble, in the case N = 100, and spectral 
endpoints ±2, at fj, = 1. 

FIG. 5. The double line diagrams for terms in the expansion of G(z) which have four y?'s in 
them. 



FIG. 6. Wigner class versus invariant class vertices. 

FIG. 7. Diagram expansions in the limit of large iV for: (a) the one-point Green function, and 
(b) the sum of one-point irreducible graphs. 

FIG. 8. Various connected quark loops with interacting gluons. Graph (a) is of order iV°, while 
(b) and (c) are of order iV -1 , and (d) is iV -2 . 



26 



FIG. 9. Graphical representation of the expansion of the "undressed" two-point connected 
Green function 

FIG. 10. (a) Diagrams representing the undressed Gaussian 3-point connected Green function 
for an even ensemble distribution, and (b) additional diagrams arising when odd moments are 
allowed. 



27 



